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APPROXIMATIONS FOR NEUTRON EMISSION SPECTRA FROM 
PROTON COLLISIONS BETWEEN 20 AND 500 MeV 
ON NUCLEI OF A > 12 


INTRODUCTION 


The approximations given in this report are based on the intranuclear-cascade data 
generated by H.W. Bertini of Oak Ridge National Laboratory (ORNL) [1], R.G. 
Alsmiller, Jr., and associates have made detailed fits of this data as described in 
References 2 and 3. The present effort was undertaken for those space shielding 
applications where the neutron spectrum may be chosen less accurately but may be more 
easily integrated into a larger computer transport code with access to a limited computer 
memory. In addition, the present fits are more satisfying because the number of 
coefficients for the data has been reduced by a factor of more than 30 ovet those of 
Reference 2, while also providing an analytical expression of the neutron emission spectra 
as a function of neutron energy, proton energy, and target A number. In Reference 3, 
Alsmiller required an elaborate double interpolation formula between ten elements and 
nine energies in order to find a neutron energy spectra formula for a given proton energy 
and target, whereas the approximate analytical expression given here completely defines 
the neutron spectra for any given proton energy and target in the range of validity. 

Of course, one should suspect that if the number of coefficients is reduced by a 
factor of 30, then the accuracy of the fits might vary accordingly. However, this was not 
the case. The difference between Alsmiller’s fits and the present work over a neutron 
spectrum was usually quite small, with a few bad exceptions, notably evaporation 
neutrons from carbon and oxygen. The agreement is better than was observed between 
the Bertini intranuclear cascade calculations and experiment. In order for the reader to 
evaluate the accuracy of the present work, as compared to Alsmiller [2,3], graphs are 
given at the end of this report for all energies and target nuclei generated by Bertini and 
given in References 2 and 3. 


USE OF CURRENT FITS 


This report provides analytical fits for two specific neutron source groups. The 
first group is the proton -induced evaporation neutrons which are integrated over all 
directions (0 to 180 deg). By dividing the analytically fit neutron energy spectrum 
by 4ir , one can derive a very reasonable estimate of the true isotropic emission spectra 
in units of neutrons/(MeV-steradian). The second group of neutron spectra is the 
proton-induced cascade neutrons which are also integrated over all directions (0 to 180 
deg). However, since cascade neutrons are fairly directional, one must use judgment in 



deriving a reasonable angular emission spectra from the given fits. If one needs an 
accurate angular emission spectra, then Alsmiller’s [2,3] data are recommended since he 
has given fits for four angular intervals (0 to 30 deg, 30 to 60 deg, 60 to 90 deg, 90 to 
180 deg). However, many good approximations for thin shields can be made using the (0 
to 180 deg) fits given here, if one can assume a straight-ahead approximation (i.e., assume 
all cascade neutrons are emitted in the same direction as the incident proton). This 
should provide a conservative estimate for normal incident protons on a slab or isotropic 
incident protons on the surface of a spherical shell with detectors at the center of the 
sphere. 

No detailed attempt will be made to explain how the data should be used except 
for the following example for neutron flux estimation: 


0(Ejsj , x) 


f E (Ep , A) 0^(Ep) F(Ej^ I Ep , A , Z) dEp 

Ep in 


neutrons 

gm-MeV 


, ( 1 ) 


where 


0(Ej^ , x) = neutron emission energy spectrum created at depth x (gm/cm^ ) , 


Z 


J (Ep , A) = inelastic cross section for proton (of energy Ep on elements of 
mass number A) in units of cm^ /gm. 


2.4 X 10~" 
aO.27 


(cm^ /gm). 


0j^ (Ep) = proton spectra at depth x(gm/cm^ ) in units of proton/(cm’ - MeV) 


and 

« 


F(EjsjjEp, A , Z) = neutron emission spectra function fit in units of neutrons/ 
MeV for a proton of energy Ep having an inelastic 

collision on a target defined by mass number A and atomic 
number Z. 


2 



Observing the function 0(Ej^j,x) in equation (1), one sees that to find the 

neutron energy spectrum at depth X* , the function must be integrated over x for all 
values from 0 to X* ; thus, for thin shields 



4>(E-^ , x) T (Efs^ , X* - x) dx 


neutrons 
cm^ -MeV 


( 2 ) 


where T(Ej^j, X* - x) is a dimensionless transport kernel for the attenuation of neutron 
flux from X to X*. When X* is sufficiently large (perhaps a mean free path), a neutron 
transport code should be used. For space vehicle shielding, it is worth noting that a 
1-MeV neutron has a mean free path in aluminum of about 15 gm/cm^, a 20-MeV 
neutron has a mean free path of about 25 gm/cm^ , and a 100-MeV neutron has a mean 
free path of about 45 gm/cm^ . If the shield contains hydrogen, the value of X* should be 
fairly small, perhaps less than 0.5 mean free path. Methods for performing the integration of 
equation (2), and especially the choice of a transport kernel will be left to the reader. 


ANALYTICAL APPROXIMATIONS FOR NEUTRON EMISSION SPECTRA 


The following pages give the neutron emission spectra for an incident proton of 
energy Ep on targets of A > 1 2. Table 1 gives the derived cascade neutron spectra 

formula for the angular interval of 0 to 180 deg. However, most of the cascade neutrons 
are emitted in the forward direction (0 to 60 deg) as measured from the proton 
direction. 

Table 2 gives the derived evaporation neutron spectra formula for the angular 
interval of 0 to 180 deg. Here it can be assumed that the neutrons are emitted 
isotropically from the excited nucleus. 

The original data generated by Bertini [1] and fit by Alsmiller, et al. [2,3], 
consisted of ten target nuclei and nine incident proton energies from 25 to 400 MeV. 
Since this report attempts to simplify the original ORNL curve fits, it seems appropriate 
to show a comparison of our results for each energy and target element that was 
originally fit by Alsmiller, even though the number was extremely large (90). Also, many 
of the neutron spectra shown here are not depicted in other reports, even though ORNL 
has extensively published most of their results. 

Figures 1 through 1 0 depict the cascade neutron spectra for ten target nuclei with 
the nine proton energies shown per figure. It should be noted that the spectral data are 
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TABLE 1 . CASCADE NEUTRON ENERGY SPECTRUM FOR PROTONS 
INCIDENT ON NUCLEI OF A> 12 

F(En|Ep , a , Z) = [Bi + B2(A-Z) + B3(A-Z)2] exp(B4t + Bjt" + B^t^); 
where 

t= ; e < Exj < E„ ; 20 < E„ < SOOMeV , 

p _ ’ rs p ’ P 

p 

e = -3.018(-l) + 4.265(-l) +5.591(-3) * , 

A 1/3 A^/3 

e = minimum neutron energy, 

Ej^ = cascade neutron energy = e + t(Ep - e) ; 0 < t < 1 , 

Ep = incident proton energy, 

A,Z = mass number and atomic number of target nuclei, 

Bi = + 5.763(-2) - 1.795(-4)Ep + 2.194(-7)Ep^ , 

Ba = - 9.499(-4) + 1.805(-5)Ep + 9.225(-8)Ep^ + 2.582(-10)Ep^ - 2.986(-13)Ep'* , 

Bj = 4.513C-6) - 1.0603(-7)Ep + 6.406(-10)Ep2 - 2.022(-12)Ep^ + 2.481(-15)Ep^ 

B 4 = - 2.272 - 1.2283(-l)Ep + 5.288(^)Ep^ - 1.1356(-6)Ep^ + 9.691(-10)Ep^ , 

Bs = + 7.845 + 1.040(-l)Ep - 2.260(-4)Ep^ + 1 ,1429(-7)Ep3 , and 

Be = - 8.007 - 6.950(-3)Ep - 7.495(-5)Ep^ + 2.0826(-7)Ep^ . 

* a.bc(Y) = a. be x 10^ 
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TABLE 2. EVAPORATION NEUTRON ENERGY SPECTRUM FOR PROTONS 
INCIDENT ON NUCLEI OF A > 12 


C(Ef,|Ep,A,Z) = C|iexp[C2;(^)*C3,(^j ] ; 


■ Neut. 
MeV 


where 


A = ^ ; 0 < Ex, < A ; 20 < < 500 MeV ; i = 1 , 2 

C4I * CsjEp N P 


Ej^ = evaporation neutron energy, 

Ep = - incident proton energy, 

A = maximum neutron energy, 

A,Z = mass number and atomic number of target nuclei. 


When i = 1 


A > 65 


Cl 1 = -t- 2.46 X 10'' + 5.5 X 10'^ (A-Z) + 5.2 x 10'* (A-Z)^ 

C 21 = - 1.87+ 1.12 X 10'^' A 

C 31 = + 5.305 x lO'* -9.646x 10"' (A-Z)'/’ 

C 41 = -t 1.92 + 3.0 X 10"’ (A-Z) 


Cs 1 = -t 3.2766 X 10"’ -t 9.874 x 10"’ (A-Z) 


'/3 


When i = 2 12 < A < 65 

Cl 2 = -1-7.0 X 10"’ (A-Z)-t 1.8 X 10"^ (A-Z)’ 

C 22 = - 7.3594 -t 1.882 (A-Z)'/’ 

C 32 = -t 6.479- 3.774 (A-Z)'/" 

A’ 

- 93.0 -t- 0.3 16 A’ 


Cs 2 = 4.39 X 10"’ + 3.40 x 10"^ A 






plotted as a function of t in the interval of 0 to 1 . This permits a simple comparison of 
all data on the same scale. The t scale is converted to the Ejj scale by the relationship 

Ej,^ = e + t (Ep - e) , where the variables are defined in Table 1 . 

Figures 1 1 through 1 9 present the evaporation neutron spectra plotted as a 
function of the neutron energy. Each of the nine figures is for a single incident proton 
energy (25 MeV-400 MeV) on each of eight target elements from oxygen through lead. 

Figures 20 through 22 show the evaporation neutron spectra for protons on 
carbon and uranium for the nine incident proton energies. The evaporation neutron 
spectra for 25-MeV protons on carbon (Fig. 20) was the poorest fit of all sets. The 
second poorest is 25-MeV protons on oxygen (Fig. 1 1). 


. 6 





(A3|M / SNOdinSN) d 
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Figure 1 . Cascade neutron spectra from carbon, A = 1 2, Z 
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Figure 2. Cascade neutron spectra from oxygen, A - 16, Z 







(A9N/SNOUin3N) d 



9 


Figure 3. Cascade neutron spectra from aluminum, A = 27, Z = 13. 
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Figure 4. Cascade neutron spectra from chromium, A = 52, Z = 24. 
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Figure 5. Cascade neutron spectra from copper, A = 65, Z = 29. 
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Figure 7. Cascade neutron spectra from cerium, A = 140, Z = 58. 
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Figure 8. Cascade neutron spectra from tungsten, A = 184, Z = 74. 
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Figure 9. Cascade neutron spectra from lead, A = 207, Z = 82. 
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Figure 10. Cascade neutron spectra from uranium, A = 238, Z = 92. 





Figure 1 1 . Evaporation neutron spectra for 25-MeV protons incident on nuclei 
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Figure 12. Evaporation neutron spectra from 50-MeV protons incident on nuclei. 
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Figure 13. Evaporation neutron spectra from 100-MeV protons incident on nuclei 
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Figure 14. Evaporation neutron spectra for 150-MeV protons incident on nuclei. 
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Figure 15. Eva] 
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Figure 16. Evaporation neutron spectra for 250-MeV protons incident on nuclei. 
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Figure 18. Evaporation neutron spectra for 350-MeV protons incident on nuclei. 
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Figure 
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Figure 20. Evaporation neutron sped 
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Figure 21 . Evaporation neutron spectra from carbon, A = 1 2, Z = 6, and uranium, A = 238, Z - 92. 
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Figure 22. Evaporation neutron spectra from carbon, A = 12, Z = 6, and uranium, A = 238, Z = 92. 
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